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Abstract. Polymer translocation in three dimensions out of planar confinements 
is studied in this paper. Three membranes are located ed z = —h, z = and 
z = hi. These membranes are impenetrable, except for the middle one at z = 0, 
which has a narrow pore. A polymer with length N is initially sandwiched between 
the membranes placed at z = —h and z — and translocates through this pore. We 
consider strong confinement (small h), where the polymer is essentially reduced to a 
two-dimensional polymer, with a radius of gyration scaling as R''g^^ ~ N'^^^\ here, 
i/2D = 0.75 is the Flory exponent in two dimensions. The polymer performs Rouse 
dynamics. Based on theoretical analysis and high-precision simulation data, we show 
that in the unbiased case h = hi, the dwell-time scales as N'^+'^^d ^ [n perfect 
agreement with our previously published theoretical framework. For hi = oo, the 
situation is equivalent to field-driven translocation in two dimensions. We show that 
in this case scales as N'^''2d ^ in agreement with several existing numerical results in 
the literature. This result violates the earlier reported lower bound N^^'^ for for 
field-driven translocation. We argue, based on energy conservation, that the actual 
lower bound for is N'^'^ and not N^'^'^. Polymer translocation in such theoretically 
motivated geometries thus resolves some of the most fundamental issues that are the 
subjects of much heated debate in recent times. 
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1. Introduction 

Polymer translocation through narrow pores in membranes is an active field of research 
in recent times: as a cornerstone of many biological processes, and also due to its 
relevance for practical applications. Molecular transport through cell membranes is 
an essential mechanism in living organisms. Often, the molecules are too long, and 
the pores in the membranes too narrow, to allow the molecules to pass through as a 
single unit. In such circumstances, the molecules have to deform themselves in order to 
squeeze — i.e., translocate — themselves through the pores. DNA, RNA and proteins 
are such naturally occurring long molecules [1-5] in a variety of biological processes. 
Translocation is also used in gene therapy [6,7], in delivery of drug molecules to their 
activation sites [8], and as a potentially cheaper alternative for single-molecule DNA or 
RNA sequencing [9,10]. 

In theoretical studies of translocation, the membrane is usually a stationary object 
that does not show any movement or fiuctuations, and this also holds for the pore in it, 
through which translocation occurs. The polymer is usually simplified to a sequentially 
connected string of monomers. A central quantity in these theoretical studies is the 
so-called dwell time r^, which is the time the pore remains blocked during a translocation 
event (Fig. [1]). 



The early theories of translocation were constructed in the spirit of mean-field [14], 
and that too for phantom polymers, wherein translocation is quantified by a Fokker- 
Planck equation for first-passage over an entropic barrier in terms of a single "reaction 
coordinate" s. Here s is the number of the monomer threaded at the pore (s = 1, . . . , A^), 
see Fig. [H These mean-field type theories apply under the assumption that every 
translocation step is slower than the equilibration time-scale of the entire polymer. 
Some years ago, this assumption was questioned [15,16], wherein the authors provided 
lower bounds for for three generic situations for phantom as well as self-avoiding 
polymers in the absence of hydrodynamical interactions: 

(a) Unbiased translocation (i.e., translocation in the absence of any driving field or 




Figure 1. Pictorial representation of a translocation event, with the polymer shown 
before, during and after translocation. We number the monomers starting with the 
end monomer on the side it moves to. The number of the monomer located within the 



pore IS s. 
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force on the polymer), for which > TRouse, with XRouse ~ N^+'^'^ being the Rouse 
time, the longest time-scale in the dynamics of the polymer; 

(b) Translocation driven by a field E, acting on the polymer only at the pore, for which 
it was shown that > N^'^^/E; 

(c) Translocation effected by a pulling force F at the head of the polymer, for which 
it was shown that > N'^/F. 

Here, u is the Flory exponent: in three dimensions u = Z/3D — 0.588 and in two 
dimensions u = z/2d = 0.75. Accompanying numerical studies led the authors to 
also suggest that the lower bounds indeed provide the correct scalings for r^; and 
based on these results, they concluded that the dynamics of translocation, in (a-c), 
is anomalous [15,16]. 

Subsequent numerical studies, however, did not immediately settle the scaling for 
Td with A^. In Tables [1] and [2] we present a summary of results on the exponent for the 
scaling of with (all results quoted are for self-avoiding polymers in the absence of 
hydrodynamical interactions in the scaling limit). 

(i) Unbiased translocation: 



authors 


two- dimensions 


three-dimensions 


Chuang et al. [15] 


2.5 




Luo et al. [17] 


2.50 ±0.01 




Wei et al. [18] 


2.51 ±0.03 


2.2 


Klein Wolterink et al. [19] 




2.40 ±0.05 


Milchev et al. [20] 




2.23 ±0.03 


Dubbeldam et al. [21] 




2.52 ±0.04 


Panja et al. [26] 




2±Z/3D 


This work 


2 ± Z/2D 





Table 1. Existing results on the exponent for the scahng of Td with N for unbiased 
translocation. 
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(ii) Field-driven translocation: 



authors 


two-dimensions 


three-dimensions 


Kantor et al. [16] 


1.53 ±0.01 




Luo et al. [22] 


1.72 ±0.06 




Cacciuto et al. [23]^ 


1.55 ±0.04 




Wei et al. [18] 




1.27 


Milchev et al. [20] 




1.65 ±0.08 


Dubbeldam et al. [24] 




1.5 


This work 


lower bound 2i/2d; 
observed 1.5 





Table 2. Existing numerical results on the exponent for the scaling of with N for 
field-driven translocation. The superscript f indicates that we discuss this paper in 
detail in the following paragraphs. 

At a theoretical level, the lack of consensus on the scalings of can easily be 
attributed to the fact that none of the works in (i) and (ii) relates the numerically 
observed scaling of to the well-known dynamical features of polymers in a satisfactory 
manner. Wc do note here that Refs. [21,24,25] proposed to link the observed scalings of 
Td to anomalous dynamics of translocation via a fractional Fokker-Planck equation; but 
how this equation can be derived from the microscopic dynamics of a single polymer, as 
well as the assumptions underlying the form of this equation, remain entirely unclear. 

In the recent past, this lack of consensus prompted us to investigate the microscopic 
origin of the anomalous dynamics for unbiased translocation. We set up a theoretical 
formalism, based on the microscopic dynamics of the polymer, and showed that the 
anomalous dynamics of translocation stem from the polymer's memory effects [26,27], in 
the following manner. Translocation proceeds via the exchange of monomers through the 
pore: imagine a situation when a monomer from the left of the membrane translocates 
to the right. This process increases the monomer density in the right neighbourhood of 
the pore, and simultaneously reduces the monomer density in the left neighbourhood of 
the pore. The local enhancement in the monomer density on the right of the pore takes 
a finite time to dissipate away from the membrane along the backbone of the polymer 
(similarly for replenishing monomer density on the left neighbourhood of the pore). The 
imbalance in the monomer densities between the two local neighbourhoods of the pore 
during this time implies that there is an enhanced chance of the translocated monomer 
to return to the left of the membrane, thereby giving rise to memory effects. The 
ensuing analysis enabled us to provide a proper microscopic theoretical basis for the 
anomalous dynamics, leading us to conclude that scales as ~ N^+^sd for unbiased 
polymer translocation in three dimensions [26,27]. 

In Refs. [26, 27] we also showed that the theory presented in Ref. [21] is not 
correct (which casts serious doubts about the correctness of a related theory presented 
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in the theoretically related paper Ref. [24]), but for unbiased translocation in three 
dimensions, the numerical result ~ jY2.52±o.o4 ^21], obtained by the use of a polymer 
model very different from ours, is consistent with Td scaling as ~ N'^+^^° . Two of us 
subsequently extended the theoretical formalism Refs. [26, 27] to analyze translocation 
by pulling the head of the polymer by a force F, leading to the theoretical derivation 



The purpose of this paper is to push the theoretical formalism of Refs. [26-28] 
further to study translocation in three dimensions out of planar confinements for 
polymers performing Rouse dynamics. Clearly, confinement reduces the number of 
configurational states available to the polymer, reducing the polymer's entropy and 
thereby increasing the polymer's free energy [29,30]. If the polymer is allowed to escape 
from the confinement through a pore, then it will translocate out of the confinement, 
and the free energy difference between the confined and the free state of the polymer 
will drive translocation. While confinement plays an important role for polymers in 
various biological processes [31], our interest in translocation out of planar confinement 
in this paper stems more from a theoretical point of view — we aim to demonstrate that 
our theoretical formalism [26, 27] works beautifully also in two dimensions. We divide 
the three-dimensional space into two parts: 2; > and ^ < by a membrane placed 
at 2 = 0. This membrane is impenetrable to the polymer except for a narrow pore. 
We then place two more parallel completely impenetrable membranes at 2; = —h and 
z = hi. The polymer is initially sandwiched between the membranes placed at z = —h 
and z = 0. We only consider strong confinement of the polymer, i.e., h -C Rf°\ with 
the radius of gyration for the polymer scaling in the present case as ~ N^^° . We 
study translocation out of planar confinement for two separate cases: (1) hi = h, and 
(2) hi = 00. Our system for /ii = cxd is shown below in Fig. [2j 



for Td ~ N^F [28]. 



monomer 1 




monomer 



Figure 2. Our system, and a snapshot of a translocating polymer out of planar 
confinement, with hi = 00. We use strong confinement, for which h Rg^^\ where 
Rf^^ is the radius of gyration for the polymer, scaling as ~ TV^sd _ 



We substantiate our theoretical analysis with extensive Monte Carlo simulations, 
in which the polymer performs single- monomer moves. The definition of time is such 
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that single-monomer moves along the polymer's contour are attempted at a fixed rate 
of unity, while moves that change the polymer's contour are attempted ten times less 
often. Details of our self-avoiding polymer model in 3D can be found in Refs. [27,32]. 

At strong confinements, i.e., with h <^ N'^^^ the confined segment of the polymer 
essentially behaves as a two-dimensional polymer. We demonstrate this below for 
h = 3 (lattice units). We tether one end of the polymer of length N/2 at the pore 
on the membrane at 2; = 0, confine the polymer between the plates at z = and 
z = —h, and measure the end-to-end distance Re of the polymer in equilibrium, as 
well as the equilibrium correlation function of the end-to-end vector. That scales 
as {N/2Y^'° and the equilibrium correlation function of the end-to-end vector behaves 
as exp[— t/rRousc(2D)] are demonstrated below in Tableland Fig. [3] respectively, with 
TRouse(2D) being the Rouse time in two dimensions, scaling, for a polymer of length N/2, 
as (Ar/2)i+2^2D_ 



N/2 


(Re) 


(i?e)/(iV/2)^- 


100 


10.248 


0.324 


150 


13.759 


0.321 


200 


17.135 


0.322 


250 


20.242 


0.322 


300 


23.171 


0.321 


350 


26.213 


0.324 


400 


28.872 


0.323 



Table 3. Average end-to-end distance for a polymer, confined between the membranes 
at z = and z = —h, with one end tethered at the pore on the membrane at z = 0, 
for h = 3. The angular brackets denote an average over 10,000 runs for each N. 

The confinement of a polymer of length N between the two planes at 2; = and 
z = —h is accompanied by an entropic (or free energ y) cost of AF ~ iV/i-i/'^3D [30]. 
Thus, in the case that hi > h, the initial state of the polymer is entropically less 
favourable, and the polymer will escape through the pore to the wider space between the 
membranes at z = and z = hi. This process is analogous to field-driven translocation, 
with a field strength ~ [/i^^/'^-'d — h-^^^^'^''^^]. Although this analogue has been correctly 
identified in Ref. [23], the interpretation of therein is not correct: Ref. [23] used z/gn to 
describe the size of the confined polymer, while Table [3] and Fig. [3] clearly show that the 
scaling of the size of the confined polymer with length is characterised by the exponent 

The proper interpretation of the numerical result of Ref. [23], therefore, clearly 
violates the lower bound A^^+^2d^ provided by Ref. [16] for r^. Table [2] shows more 
independent numerical evidence that the scaling of for field-driven translocation falls 
far short of A^i+'^ao. This raises serious doubts about the theoretical lower bound A^^"^'" 
for Trf for field-driven translocation argued in Ref. [16]. 
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Our main results in this paper are two-fold. First, for hi = h, the entropic drive 
is absent; the translocation dynamics reduces to that of an unbiased translocation for a 
two-dimensional polymer. Based on theoretical analysis and high-precision simulation 
data, we show that the dwell time for a polymer of length scales as ~ N'^'^'^^o 
for hi = h, in perfect agreement with our previous results [26,27]. In this paper, we 
actually go one step further than Refs. [26,27] to show that the probability distribution 
of the dwell time, P{Td) for hi = h has a scaling form P{Td) ~ V{Td/N'^^''^°) / N^+''''° . 
Secondly, for field-driven translocation, we argue, based on conservation of energy, that 
the lower bound for for field-driven translocation is given by N'^'^ in the absence of 
hydrodynamics. Using the analogue between translocation out of planar confinement for 
hi = oo and field-driven translocation in two dimensions, we demonstrate numerically 
that Td ~ A^^'^aD^ with our numerical results being consistent with those of Refs. [16] 
and [23]. Study of polymer translocation in these (theoretically motivated) geometries, 
therefore, is a fine test case for the fundamental physics governing translocation 
dynamics. 

1.1 



0.7 

m 

0.3 

0.0 
-0.1 

2.4e5 4.8e5 7.2e5 

(300)'"''f/A^'' 

Figure 3. End-to-end equilibrium correlation function c{t) for a polymer with one 
end tethered at the pore on the membrane at z = 0, for h = 3. The data for each N 
are averaged over 256 realisations. 

This paper is organised as follows. In Sec. [2] we discuss a method to measure 
^{t), the component of the polymer chain tension at the pore, perpendicular to the 
membrane. In Sec. E] we analyze the memory effects in the imbalance of the 

polymer chain tension across the pore. In Sec. H] we discuss the consequence of these 
memory effects on unbiased translocation, i.e., for the case hi = h. In Sec. EJwe derive 
the lower bound N"^^ for for field-driven translocation, and discuss the consequence 
of this lower bound and the memory effects on translocation out of confinement for the 
case hi = oo. We finally end this paper with a discussion in Sec. El 




Polymer Translocation out of Planar Confinements 



8 



2. Chain tension at the pore perpendicular to the membrane 

A translocating polymer can be thought of as two segments of polymers threaded at the 
pore, while the segments are able to exchange monomers between them through the pore. 
In Ref. [26] we developed a theoretical method to relate the dynamics of translocation to 
the imbalance of chain tension between these two segments across the pore. The key idea 
behind this method is that the exchange of monomers across the pore responds to the 
imbalance of chain tension in its turn, (f){t) adjusts to v{t), the transport velocity of 
monomers across the pore. Here, v{t) = s{t) is the rate of exchange of monomers from 
one side to the other, where [s{t) — s(0)] is the total number of monomers translocated 
from one side of the pore to the other in time [0,t]. In fact, we noted that [s(t) — s(0)] 
and (j){t) are conjugate variables in the thermodynamic sense, with playing the role 
of the chemical potential difference across the pore. 



0.92 

0.72 
0(f = 0) 

0.52 
0.32 




N 1.1 
V 







. / 




I 
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d 



10 



1.1 



1.2 



1.3 



,(4) 



< Z (f = 0)> 



Figure 4. {Z^^\t ~ 0)) vs. <i>(i = 0) demonstrating the linear relationship between 
the two, for N — 100 and h — 2,3,4,5,6 and 8 respectively. The angular brackets 
for = 0)} indicate an average over 2,400 polymer realisations. The data for 

= 0) are obtained 2,400 polymer realisations as well. The solid line corresponds 
to the linear best-fit. Inset: {Z'^'^\t = 0)) as a function of h. 



By definition, (pit) = ^z>o{t) — ^z<o{t) where $2>o(^) and $2<o(t) are respectively 
the chain tension (or the chemical potential) on the 2; > and the z < side of the 
pore. Consider a separate problem, where we tether one end of a polymer to a fixed 
membrane, yet the number of monomers are allowed to spontaneously enter or leave the 
tethered end, then we have 

^|^^ = exp[<|.(t)A^T], (1) 

where Wt{ > +) [resp. Wt{-\ > —)] is the rate that a monomer enters (resp. leaves) 

the polymer chain through the tethered end at time t. Note that tethering the polymer 
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while allowing monomers to enter or leave the polymer at the tethered end is precisely 
the case that translocation represents. 

For a translocating polymer out of confinement between the membranes at ^ = 
and z = —h, note that at t = 0, it is easy to use Eq. ([1]) to measure the chain tension 
for both segments at the pore [^{t = 0) in our notation], since under these conditions, 
we also have the relation that 

P_ Wt=o{- ^+) = P+ m=o(+ -) , (2) 

where P_ (resp. P4.) is the probability that the z < (or the 2; > 0) polymer segment 
has one monomer less (resp. one extra monomer). Equations ([1]) and ([2]) together yield 
us 

^t = 0) = kBT\n^. (3) 

The chain tension as obtained from Eq. ([3]) is linearly related to the distance 
of the centre-of-mass of the first few monomers along the polymer's backbone, at the 
immediate vicinity of the pore, at least in our simulations. This is shown in Fig. IU 
where for a tethered polymer of length N = 100, the average distance {Z^^\t = 0)) 
of the centre-of-mass of the first 4 monomers along the polymer's backbone from the 
membrane, counting from the tethered end of a polymer, is plotted versus the chain 
tension $, for a variety of h values. Within the error bars, all the points in Fig. Hlfall 
on a straight line, implying that $ is very well-proxied by {Z^'^^). Since measurements 
of the chain tension via Eq. ([3]) are much more noisy than measurements of {Z^'^^), we 
will use the latter quantity as a measure for the chain tension. 

3. Memory effects in the ^-component of the chain tension at the pore 

In the case of unbiased polymer translocation, we have witnessed in Refs. [26,27] that 
the memory effects of the polymer give rise to anomalous dynamics of translocation. 
We argued [26,27] that the velocity of translocation v{t) = s{t), representing monomer 
current, responds to 0(t), the imbalance in the monomeric chemical potential across the 
pore acting as "voltage". Simultaneously, (pit) also adjusts in response to v(t). In the 
presence of memory effects, they are related to each other by (pit) = (pt=o + Jq dt'ii{t — 
t')v(t') via the memory kernel /i(t), which can be thought of as the (time-dependent) 
"impedance" of the system. 

In this section, following Refs. [26, 27] we determine the memory kernel n{t) to 
describe the dynamics of translocation out of (strong) planar confinements. Note that 
for the case hi = h, there is an obvious symmetry between the polymer segments 
confined within the parallel plates below and above the z = plane, implying that the 
corresponding memory kernels denoted by /ii>o(^) and /xi<o(^) are the same. For the 
case hi = 00, we already know the form of /ii^o(^) from Ref. [26,27]. In fact, in Ref. [26] 
we determined /ii'^o(^) ^y injecting p monomers into the tethered end of an equilibrated, 
tethered polymer of length N/2 — p (bringing the total length to N/2), and proxying 
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(pit) by the average distance of the centre-of-mass of the first 4 monomers {Z^^\t)) from 
the membrane. We found 

/^i>o(^) ~ t '^'^^^ exp[-t/rRousc{3D)] , (4) 

with rRouse(3D) is the Rouse time for a polymer of length N/2, i.e., rRouse(3D) ~ 
(Ar/2)i+2-«D. 

Following Refs. [26,27], here we compute /ii<o(^)! the memory effect of a polymer 
of length N/2 with one end tethered to the pore in the membrane placed at z = 0. 
Indeed, the expression for /i^'^^Q(t), as we derive below, is given by 

/^2<o(^) ~ t '^'^'^ exp[-t/rRouse(2D)] • (5) 



10" N ^ 



: is 




I A I I 

10 10 10 10 10 10 
t 

Figure 5. (Z('')(oo) - Z^'^^t)) for iV/2 = 250, for /i = 3 (pluses) and /i = 2 (squares). 
The solid line corresponds to the power law i^(i+'^2D)/(i+2i'2D), Note that the data for 
h = 2 obeys a cleaner power law t^(^+'^2D)/(i+2i'2D) than the data for h — 3: this we 
can expect, as the polymer resembles a two-dimensional polymer more for h — 2 than 
for h = 3. The angular brackets correspond to 3, 200, 000 polymer realisations. 



While the Rouse relaxation exp[— t/rRouse(2D)] can be easily justified based on Fig. [HI 
the value of a for /ii'^o(^) ~ exp[— t/rRouse(2D)] is obtained by following the procedure 
of Refs. [26,27]. The value of a depends on the relaxation properties following the event 
of injecting, say, p extra monomers at the tether end, just like extra monomers add to 
(or get taken out of) the polymer segment confined within z < during translocation. 
Given the exp[— t/rRouse(2D)] behaviour of Fig. [3l we anticipate that by time t after the 
extra monomers are injected at the tethered point, the extra monomers will come to a 
steady state across the inner part of the polymer up to Ut ~ ti/(i+2'^2D) monomers from 
the tethered point, but not significantly further. This internally equilibrated section of 
Ut+p monomers extends only r{nt) ~ n^^° , less than its equilibrated value [uf + pY^^ , 
because the larger scale conformation has yet to adjust: the corresponding compressive 
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force from these Ut+p monomers is expected by standard polymer scaling [30] to follow 
f /{kBT) ~ 5r{nt)/r'^{nt) ~ i^2dP/ [ntr{nt)] ~ t-(i+''2D)/(i+2!^2D)^ for p <^ As was the 
case in Refs. [26-28] , we expect that the chain tension at the pore behaves linearly with 
the force /, leading to a = (1 + z/2d)/(1 + 2z/2d) = 0.7. 

We have confirmed this picture by measuring the impedance response through 
simulations. In Fig. |H we have shown that the centre-of-mass of the first few monomers 
is an excellent proxy for chain tension at the pore and we assume here that this further 
serves as a proxy for 5$. Based on this idea, we track by measuring {Z^^\t))^ 

in response to the injection of extra monomers near the pore at time t = 0. Specifically 
we consider the equilibrated segment of the polymer confined within 2; < 0, of length 
N/2 — 5 (with one end tethered at the pore), adding 5 extra monomers at the tethered 
end of the polymer segment at time t = 0, corresponding to p = 5, bringing its length 
up to N/2. Using the proxy {Z^^\t)), we then track (S^^'^^^^t)) . At short times, these 
p monomers quickly expand in the 2;-direction (e.g., t < 400 in Fig. [5]). Only after they 
"feel" the impenetrable membrane ai z = —h, they turn around to expand along the 
xy-plane: indeed, the agreement between the latter and the theoretical prediction of 
a = (1 + 1^20)/ (1 + 2z/2d) for t > 400, for N/2 = 250, and for /i = 2 and 3, can be seen in 
Fig. [5l Note that the sharp deviation of the data from the power law at long 

times is due to the asymptotic exponential decay as exp[— t/rRousc(2D)] of {S^^^^'^\t)) at 
long times. See also Fig. 2 of Ref. [26] in this context. 

4. The case hi = h: consequence of the polymer's memory effects on 
unbiased translocation 

For the case hi = h the condition fi^^yo{t) = /ii<o(^) = f^'"'^\t) for the polymer segments 
above and below the plane z = allows us to write (j){t) = (f)t=o 

as presented in the starting paragraph of Sec. [31 Supposing a zero-current equilibrium 
condition at time t = 0, this relation can be inverted to obtain v{t) = dt'a{t — t')(f){t'), 
where a{t) can be thought of as the "admittance". In the Laplace transform language, 
jl'^^\k) = d~^{k), where k is the Laplace variable representing inverse time. Via 
the fluctuation-dissipation theorem, they are related to the respective autocorrelation 
functions as //(^^(t - t') = {(f){t)(f){t'))^=o and a{t - t') = {v{t)v{t'))^=o- 

As explained in the introduction, translocation for hi = h is unbiased, for 
which, having shown that jj.'^^^t) ~ t 1+2^20 exp[— t/rRousc(2D)]5 we expect [26, 27] 
that the translocation dynamics is anomalous for t < trousc{2D), in the sense that 
the mean-square displacement of the monomers through the pore, (As^(t)) ~ t^ 
for some (3 < 1 and time t < TRouse(2D), whilst beyond the Rouse time it becomes 
simply diffusive. Strictly speaking, rRouse(2D) in this expression should be replaced 
by the characteristic equilibration time of a tethered polymer with length of 0{N); 
since both scale as N^+^'^^o^ ^gg rRouse{2D) here, favouring notational simplicity. 
The value P = a = i^2j^°^ follows trivially by expressing (As^(t)) in terms of 
(translocative) velocity correlations {v{t)v(t')), which (by the Fluctuation Dissipation 
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theorem) are given in terms of the time dependent admittance a{t — t'), and hence 
inversely in terms of the corresponding impedance. In other words, up to the 
Rouse time, the squared displacement as a function of time is subdiffusive, following 
(As^(t)) ~ ti+2''2D. Consequently, at the Rouse time rRouse(2D) ~ A^i+2!/2d^ ^\^q squared 
displacement scales as (As^[rRouse(2D)]) ~ N^+'^^o_ Beyond the Rouse time, there are 
no memory effects and the squared displacement increases linearly in time: (As^(t)) ~ 
{(As^[rRouse(2D)])/TRousc{2D)}^ ~ N~'^^°t. Based on the criterion for unthreading, i.e., 
unthreading occurs when ^{As^{Td)) ~ N, one then obtains Td ~ iV^+'^aD [26,27]. 



N 








30 


2439 


0.2114 


0.4948 


40 


5176 


0.2034 


0.5115 


50 


9499 


0.2021 


0.5373 


60 


15684 


0.2021 


0.5624 


70 


24532 


0.2069 


0.5984 


80 


34556 


0.2018 


0.6037 


90 


47974 


0.2027 


0.6243 


100 


64755 


0.2048 


0.6476 


200 


415767 


0.1954 


0.7350 


300 


1268463 


0.1955 


0.8137 


400 


2765246 


0.1932 


0.8641 


500 


4961331 


0.1877 


0.8875 


600 


8228721 


0.1885 


0.9332 


700 


12648891 


0.1897 


0.9758 


800 


17975330 


0.1867 


0.9930 



Table 4. Median values of based on 8, 192 unthreading realisations for each N. 

For computer simulations of unbiased translocation in two dimensions, Luo et 
al. [17] reported a scaling of ~ ]yi+2i/2d. note that this exponent and our theoretical 
expectation is 10% different. To distinguish these two different exponents, we performed 
high-precision simulations to obtain the unthreading time for a number of iV-values: 
N = 30,40,50,60,70,80,90,100,200,300,400,500,600,700,800, for 8,192 realisations 
for each value of A^, with a membrane spacing of h = 3 (Table Hj). The unthreading time 
Tu in Table HI) is defined as the time for the polymer to leave the pore with s{t = 0) = N/2 
and the two polymer segments equilibrated at t = 0. Both and scale the same way, 
since Tu < < 2r„ [27]. The data of Table HI along with the effective exponent for r„ 
vs. as a power law are further shown in Fig. [6l To be able to use the full potential 
of the statistics of 8, 192 realisations for each A^, we also plot the sorted t^/N'^^^'^^ and 
Tu/A^^+^'^^D yg_ ^j^g normalised rank of the sorted values in Fig. [71 The data collapse is 
a further test that rules out the scaling as A^^+^'^ac^ reported in Ref. [15,17,18]: in 
fact the left panel of Fig. [7] suggests that P{Td), the probability distribution for the dwell 



Polymer Translocation out of Planar Confinements 13 




Figure 6. Scaling of t„ with N: r„ data of Table |4] are represented by the black line 
with points, solid line corresponds to the scaling t„ ~ TV^+^sd _ xhe effective exponents 
(9(lnr„)/a(ln A^) for N = 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700 are shown 
in the inset, clearly ruling out A^i+2'^2d scaling of t„ in favour of N'^+'^2d ^ 




0.5 1 0.5 

Normalized rank after sorting on x 



Figure 7. Sorted values Tu/N'^^'^'^'o and Tu / N^^'^^'^'o vs. their normalised rank, for 
8, 192 realisations per value of N . Data shown (from bottom to top in the right panel): 
N = 30, 40, 50, 60, 70, 80, 90, 100, 200, 300, 400, 500, 600, 700, 800. The left panel shows 
a much better collapse than the right panel, ruling out the Td scaling as N^+'^'^^o ^ 
as observed by Luo et al. [17]. In fact, the data collapse in the left panel suggests 
that P(T(j), the probability distribution for the dwell time has a scaling form 
P{Td) ~ 7'(Td/iV2+>^2D)/iv2+''2D^ with a scaling function V{x). 



time Td has a scaling form P{Td) ~ V{Td/ N'^'^^^^) / N'^^'^'^^ ^ with a scahng function V{x). 
These results together clearly demonstrate that the dwell time scales as Td ~ N'^+'^^'°. 
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5. for hi = oo, or equivalently, r^^ for field-driven translocation in two 
dimensions 

5.1. Lower bound for for field- driven translocation 

We now turn to the case hi = oo, which is equivalent to field-driven translocation in two 
dimensions as we discussed in the introduction. For our first stop, we notice in Table 
[2] that quite a few reported numerical results violate the lower bound A^^^'' for the 
scahng exponent of r^, suggested in Ref. [16], including that of the authors of Ref. [16] 
themselves. In light of this, below we first discuss the lower bound for field-driven 
translocation. 

The crux of the derivation of the lower bound for in Ref. [16] is that, with or 
without an applied field, the mobility of a polymer translocating through a narrow pore 
in a membrane will not exceed that of a polymer in bulk (i.e., in the absence of the 
membrane). To obtain the mobility of a polymer in bulk, the authors assumed two 
more attributes of a polymer under a driving field: 

(i) To mimic the action of a field on a translocating polymer, the field on the polymer 
in bulk has to act on a monomer whose position along the backbone of the polymer 
changes continuously in time. As a result, there is no incentive for the polymer 
to change its shape from its bulk equilibrium shape, i.e., the polymer can still be 
described by a blob with radius of gyration ~ A^'' in the appropriate dimension. 

(ii) The polymer's velocity is ~ mE, where E is the field, and m is the mobility ~ 1/A^. 

Of these two assumptions, note that (ii) is obtained as the steady state solution of the 
equation of motion of a Rouse polymer, in bulk, with uniform velocity and vanishing 
internal forces (see e.g., Ref. [30], Eq. VI. 10). We have already witnessed in many 
occasions [15,16,21,24-28] that the dynamics of translocation through a narrow pore 
is anomalous (subdiffusive), and in Sec. [3] of this paper we have seen that there are 
strong memory effects in the polymer, to the point that the velocity of translocation 
is not constant in time. The anomalous dynamics and the memory effects are crucial 
ingredients that question the validity of the lower bound A^^"'"'^ for for field-driven 
translocation. 

A lower bound for for field-driven translocation does nevertheless exist, and it 
can be obtained from conservation of energy. Consider a translocating polymer under 
an applied field E which we can assume to be acting only at the pore: monomers 
take time to translocate through the pore. The total work done by the field in time 
Td is then given by EN. In time r^, each monomer travels a distance of ~ Rg, leading 
to an average monomer velocity Vm ~ Rg/Td. The rate of loss of energy due to viscosity 
rj of the surrounding medium per monomer is given by r]v'^. For a Rouse polymer, 
the frictional force on the entire polymer is a sum of frictional forces on individual 
monomers, leading to the total free energy loss due to the viscosity of the surrounding 
medium during the entire translocation event scaling as AF ~ Nrdrjv^ = NrjRg/Td- 
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This loss of energy must be less than or equal to the total work done by the field EN, 
which yields us the inequality > rjRg/E = rjN'^'^/E [33]. 

5.2. Td for the case hi = oo 

If we follow the procedure due to two of us in Ref. [28] to calculate for the case 
hi = oo via the memory kernels discussed in Sec. [H then we would adopt the following 
route. For the translocated part of the polymer (in the space z > 0), the memory 
kernel takes the form of Eq. (jl]), while the translocating part of the polymer (in the 
space z < and z > —h) the memory kernel takes the form of Eq. (I5l). Of these, the 
magnitude of the exponent in the power law of Eq. ([5]) is less than that of Eq. (jlj) [i.e., 
the memory effects of the translocating part of the polymer are longer-lived than that 
of the translocated part of the polymer], which implies that the relation between the 
chain tension imbalance across the pore and translocation velocity should be described 
by the equation 



via Laplace transform [28]. Furthermore, if [0(0) — 0(t')] remains a constant (not shown 
here [34]), then Eq. (CD yields v{t) ~ t-'^aD/Ci+Zi^aD)^ implying that the distance [s{t) - 

s(0)] unthreaded in time t should behave as s{t) = s(0) + / dt' v{t') ~ t(i+'^2D)/(i+2i^2D)_ 



With [s{Td) - s(0)] = N, the relation [s(t) - s(0)] ~ t^^+'^'2B)/ii+2u^B) ^^^i^^ finally yield 

The scaling Td ~ jv(i+2;^2D)/(i+i'2D) obtained through the memory kernel approach 
violates the lower bound ~ j\r2i/2D^ ^nd therefore the former cannot be the correct 
scaling for Td- Indeed a short reflection makes the issue clear. The memory kernels 
we derived in Sec. |3] are actually "static memory kernels" : the individual monomer 
velocities involved, for determining the static memory kernels (i.e., for the equilibration 
process of the polymer when p extra monomers are injected at the pore, with p -C A^), 
are small, and as a result, the energy loss due to viscosity of the surrounding medium 
is negligible. For fleld-driven translocation in two-dimensions, since the lower bound 
for Td set by the viscous energy loss overrides the expression of Td obtained from the 
static memory kernel, it is very well possible that there exists a corresponding dynamic 
memory kernel" that describes the relations between (j){t) and v{t) for a translocating 
polymer. 

We have not found a way to probe this dynamic memory kernel. We can nevertheless 
view its consequences by plotting the distance unthreaded [s(t) — s(0)] as a function 
of t. We chose h = 3 for a variety of lengths of polymers (corresponding to strong 
confinement, i.e., equivalent to translocation in two dimensions driven by a strong field) 




(6) 



leading to 
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Figure 8. The mean time (t) required to unthread a distance [s — s(0)] for h = 3, 
hi = oo, s(0) = 3iV/4 and [s — s(0)] = 5, 10, . . . , iV. The (t) values are obtained as 
average over 2, 048 reahsations for each N. The iV- values used are (from left to right): 
200, 280, 360, 400, 480, 520, 600, 680, 800. The (i)-values for N = 200 is the actual value 
of time, for the others, (t)-values for each N are separated by a factor 2 along the x- 
axis. The solid line represents a power law with exponent 1/(2z/2d), corresponding to 



for this purpose. For a polymer of length N we hold 3N/ 4 monomers between the plates 
z = and z = —h, with the other N/A monomers protruding out through the pore in 
the space z > 0, i.e., s(0) = SN/A [35]. We equilibrate both segments, and at t = we 
let translocation begin. The mean time (t) required to unthread a distance [s — s(0)] 
and the scaling of are presented in Fig. [8] and Table [5] respectively. We note that 
although A^^'^sD provides only the lower bound for r^, our simulation data, along with 
those of Refs. [16,23] suggests that A^^^^d jg indeed the correct scaling for the in two 
dimensions. 

6. Discussion 

Polymer translocation out of confined spaces plays an important role for polymers in 
various biological processes. Planar confinements, however, are mostly a theoretical 
construct. In this paper, we have studied polymer translocation in three dimensions 
out of planar confinements, and demonstrated that polymer translocation in these 
(theoretically motivated) geometries is a very interesting testcase of fundamental physics 
of translocation dynamics. 

The geometry we have considered is as follows. We have divided the three- 
dimensional space into two parts: 2; > and z < by a membrane placed at z = 0. This 
membrane is impenetrable to the polymer except for a narrow pore. We have then placed 
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N 






200 


117934 


41.6960 


280 


198497 


42.3659 


360 


294199 


43.0712 


400 


358865 


44.8581 


480 


454131 


43.1836 


520 


519663 


43.8245 


600 


653885 


44.4912 


680 


779484 


43.9586 


800 


1026083 


15.3169 



Table 5. Median values of tj, for h = 3, demonstrating tj, ~ N'^'^w , xhe values are 
based on 2, 048 translocation realisations for each N. 

two more parallel membranes at z = —h and z = hi that are completely impenetrable 
to the polymer. The polymer is initially sandwiched between the membranes placed 
at z = —h and z = 0. We have considered strong confinement for the polymer, i.e., 
h <^ Rij^°^ where R^g°^ is its radius of gyration in the bulk scaling as ~ N'''-^'° . Here, 
N is the polymer length and ~ 0.588 is the Flory exponent in three dimensions. 
Under these conditions the confined segment of the polymer essentially behaves as a 
polymer in two dimensions. If hi > h, the initial state of the polymer is entropically 
unfavourable, and the polymer escapes to the space between the membranes placed 
at z = and z = hi through the pore in the membrane placed at ^ = 0: this is 
essentially the field-driven translocation process in two dimensions. We have studied 
two separate cases: (i) hi = h, and (ii) hi = oo. For (i) the entropic drive is absent; 
the translocation dynamics reduces to that of an unbiased translocation for a two- 
dimensional polymer. Based on theoretical analysis and high-precision simulation data, 
we have shown that the dwell time r^, the time the pore remains occupied during 
translocation, scales as N'^+'^^o ^ perfect agreement with our previous results. We have 
also shown that the probabihty distribution of the dwell time, P{Td) for hi — h has a 
scahng form P(rd) ~ Vird/N^^'"'^) /N'^^''^^ . For (ii) we have shown that Td ~ N'^''^^ , 
in agreement with several existing numerical results in the literature. Here z/2D = 0.75 
is the Flory exponent in two dimensions. The result ~ jV^'^^D ^qj- (^jj^ violates 
the earlier reported lower bound 1 + A^'^ for Td for field-driven translocation. We have 
argued, based on conservation of energy, that the actual lower bound for Td is N^'^ and 
not 1 + A^''. 
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